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We propose to produce five-partite entanglement via cascaded four-wave mixing in a high-Q mi¬ 
croresonator that may become a key to future one-way quantum computation on chip. A theoretical 
model is presented for the underlying continuous-variable entanglement among the generated comb 
modes that is expansible to more complicated scenarios. We analyze the entanglement condition 
when the van Loock and Furusawa criteria are violated, and discuss the device parameters for po¬ 
tential experimental realization that may be utilized to build an integrated compact five-partite 
entanglement generator. The proposed approach exhibits great potential for future large-scale inte¬ 
grated full optical quantum computation on chip. 
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I. INTRODUCTION 

Quantum computation (QC) is expected to provide 
exponential speedup for particular mathematical prob¬ 
lems such as integer factoring [I| and quantum system 
simulation [^. However, any practical QC system must 
overcome the inevitable decoherence problem and achieve 
scalability. The traditional “circuit” QC model keeps 
quantum information in a physical system where quan¬ 
tum memory units undergo precise controlled unitary 
evolution simultaneously, leading to serious scalability is¬ 
sue. To circumvent this challenge, an “one-way” quan¬ 
tum computation model was proposed Q, where quan¬ 
tum information exists virtually in a cluster state y and 
one can perform any desired quantum algorithm by con¬ 
ducting a sequence of local measurements. With this 
approach, the most challenging part is now transferred 
from conducting the unitary operation in a large scale 
into the generation of a cluster state, or more generally, 
a universal multipartite entangled state. The aim of this 
paper is to investigate the possibility of a novel integrated 
approach for generating multipartite entangled states. 

Optical frequency combs (OFCs) have been shown to 
be capable of preserving cluster states [i,[i]. An OFC 
is a light source composed of equally spaced discrete fre¬ 
quency components, as illustrated in Fig. [T](a). Actual 
OFCs might extend to an extremely broad band with 
hundreds of frequency components [3, , each of which 

corresponds to a comb mode (marked by a mode number, 
say m). OFCs are favorable for QC for their robustness 
to decoherence since photons are less likely to inter¬ 
act with the environment compared with other physical 
systems such as atoms 

OFCs have already been utilized in many applications 
such as frequency metrology, telecommunications, opti¬ 
cal and microwave waveform synthesis, and molecular 
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spectroscopy Hi- Conventionally, OFCs are generated 
in mode-locked lasers that are usually bulky, difficult to 
^erate, and susceptible to environmental perturbations 
Q. It is recently reported that OFCs can also be gen¬ 
erated from monolithic microresonators [113 through 
cascaded four-wave mixing (FWM). 

In a high-Q microresonator with appropriate disper¬ 
sion, an intense pump wave launched into a cavity mode 
would excite four-wave mixing processes among differ¬ 
ent cavity modes via the optical Kerr effect [l^. There 
are dominantly two types of FWM, degenerate and non¬ 
degenerate, which are illustrated in Fig.[T](b). Due to the 
momentum conservation among the interacting photons, 
a degenerate process converts two identical photons in 
a same mode at m into two dissimilar photons at modes 
m—1 and m-|-I, respectively. Similarly, a non-degenerate 
process converts two photons from modes m and m -I- I 
into two new photons at modes m — I and m -|- 2. The 
iteration of these two processes thus produce an optical 
frequency comb 3, with a spectral extent determined 
by the group-velocity dispersion of the device. 

The beauty of such scheme lies in the nature of high- 
Q microresonators. First, the optical field is strongly 
confined inside a small volume, leading to significantly 
enhanced nonlinear optical interactions. Second, due to 
the exceptionally high quality factors {Q) of microres¬ 
onators, the photon life time inside the cavity is much 
longer than that in those traditional cavities so that dif¬ 
ferent frequency components have enough time to entan¬ 
gle with each other. Finally, the integrated chip-scale 
platform of microresonators exhibit great potential for 
eventually realizing a large-scale integrated full optical 
quantum computer 3- 

These facts inspire us to explore the potential of OFCs 
for producing multi-partite entangled states inside a mi¬ 
croresonator. Although two-mode quantum squeezing 
has been intensively investigated for parametric processes 
[13 - 3 1 the quantum properties of microresonator-based 
frequency comb generation has not yet been fully ad¬ 
dressed. On the other hand, there have been both 
theoretical analysis [23 - 13 experimental invest!- 
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FIG. 1. (a) The spectrum of an ideal frequency comb is dis¬ 
crete, equally spaced, and covers a wide band, (b) Energy 
level diagram of degenerate (left) and non-degenerate (right) 
FWM. 


gation [23 - [^ on photon pair generation inside mi- 
cro/nanophotonic devices. Yet all of them focused on 
the bipartite discrete-variable entanglement and the bi¬ 
partite methodology cannot be applied to the analysis 
of entanglement among three or more frequency compo¬ 
nents. In this paper, we present a theoretical model to 
describe the five-partite continuous-variable (CV) entan¬ 
glement among frequency comb modes. We solve the 
Fokker-Planck equation in P representation and then an¬ 
alyze the entanglement condition when van Loock and 
Furusawa criteria are violated. We also present an ex¬ 
periment scheme to verify our prediction. 

The rest of this paper is arranged as follows: in sec¬ 
tion |lll we propose an experiment scheme to measure the 
continuous-variable five-partite entanglement generated 
from a microresonator. From section Hill to section IVl we 
analyze how the degree of five-partite entanglement may 
vary with respect to different cavity configuration. Sec¬ 
tion IVII is for some discussion on the simulated output 
fluctuation spectra. Finally the conclusions are drawn. 


II. SYSTEM MODEL 

The experiment scheme is shown in Fig. A nar¬ 
row linewidth continuous-wave laser is used to pump the 
system. The pump beam first go through a polariza¬ 
tion controller, then passes an arrayed waveguide grating 
(AWG) which serves as a narrow bandwidth bandpass 
filter. The pump beam is coupled into the nonlinear di¬ 
electric cavity using evanescent fields. The intracavity 
photons interact through FWM and generate the OFC, 
which is later extracted through evanescent fields again. 
Once the different frequency components of OFC is sepa¬ 
rated by AWG, we will be able to analyze each component 
in Fabry-Perot (FP) analysis cavities [3l|. 

When the pump is weak, no side band is generated. As 
the pump power increases, more frequency components 
will emerge EH- In this paper, we consider the situation 
when there are all together five frequency components in 
the resonator. They are marked as s2, sl,p,il,i2 from 
lower to higher frequency in the spectrum, as plotted 
in Fig. [1] (a). In practice, the total mode number of 



FIG. 2. OFC generator with a Galcium Fluoride cavity and 
angle-polished fiber couplings. GW, continuous-wave; PC, 
polarization controller; AWG, arrayed waveguide grating. 

the comb can be controlled by engineering the device 
dispesion. In the interaction picture, the Hamiltonian of 
the system can be written into three parts, i.e.. 


H — -j- Hint F Hr^ 

Hint = 

k 


pump: Hpnmp — -j- H.C.^ 

alalukak + igh'^alalatUk 

i—i-t 


-A-inh{^ nt. n'^ 


( 1 ) 

( 2 ) 


where k,t = p, si, s2,il,i2. The interaction part of the 
Hamiltonian is derived from the usual Kerr Hamiltonian 
V = fi{g/2) : {ap + asi+aii+as 2 + ai 2 + H.c.)'^ :, where 
... stands for normal ordering, g is coupling coefficient 
which will be explained later, and e is the pump field 
that enters the resonator, which, based on non-depletion 
assumption, is described classically, 

The interaction Hamiltonian consists of three parts re¬ 
sponsible for self-phase modulation, cross-phase modu¬ 
lation, and four-wave mixing, respectively. It is easy to 
verify that the first two parts automatically vanish in the 


P representation [T^. We have also omitted some FWM 
terms indicating the physically less significant nonlinear 
process for the sake of simplicity, including: as 2 ®l 2 ®p! 
corresponding to the processes that 2p —>■ s2 + i2, is omit¬ 
ted, because the inevitable dispersion will cause farther 
separated modes be less coupled; and a]. 2 Cii 2 <^siaii, cor¬ 
responding to si + il —>■ s2 -I- i2, because the absence 
of pump indicates the process will be much weaker than 
others. All these processes might actually help the phase¬ 
locking mechanism, but are less dominant to the gen¬ 
eration of entangled states, so we will leave them for 
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future work. In the end there are only hve pairs of 
FWM terms left in the hamiltonian: a\2a1pa1i, 

al 2 aj,aji, ag 2 aliasiap, their Hermitian 

counterpart. The latter two are non-degenerate. Such 
choice of FWM processes is consistent with the previously 
mentioned assumptions on the entanglement of multipar¬ 
tite OFCs. 

The coupling coefficient g is defined as 


9 = 


712 /(WqC 

Vnl 


( 3 ) 


where n 2 is nonlinear refractive index that characterizes 
the strength of the optical nonlinearity, no is the linear 
refractive index of the material, c is the speed of light in 
the vacuum. V is the mode volume. It can be approx¬ 
imated by [s^, yet a more accurate expression 

is given in [ij]. That is to say, if we can expand the 
intracavity electric field in time domain as 


E{r,t) = iep(<)e“^''*Tp(r) -k c.c., (4) 


ically utilizing eigenmode solver software snch as WG- 
MODES [13. 

To describe an open system we introduce the loss and 
out-coupling terms, 

Lkp = -ik{‘2.akpa\-a\akp-pa\ak), (6) 

in which p is the density matrix of the hve cavity modes 
under consideration. 7 ^ = 7^9 + 7 fcc stands for the damp¬ 
ing rates for the corresponding resonator modes, which 
is the sum of two terms, ^kc is related to the ampli¬ 
tude transmission coefficients, which is experimentally 
tunable, and 7 fco is related to the intracavity absorption 
rate, determined by the material and shape of the res¬ 
onator. 


III. EQUATIONS OF MOTION FOR THE FULL 
HAMILTONIAN 

Now consider the whole physical system, where the 
nonlinear dielectric is contained inside a pumped reso¬ 
nant cavity. The master equation for the hrst hve cavity 
modes is 


where is the slow-varying envelope. Then each 

mode volume is calculated as 


V^ = 


n -1 


|Tp(r)||MU 


( 5 ) 


where T^(r) is the normalized eigenvector of the trans¬ 
verse electric (TE) or transverse magnetic (TM) modes. 

For a spherical cavity, an analytical solution for T^(r) 
is given in 0 . In our case Vq = 6.6 x 10 m^. 

For other cavities, for example, a micro-rin g re sonator, 
it might be infeasible to solve analytically [^. How¬ 
ever, it might still be possible to calculate EqjS] numer- 


dt 


i ^ 

^ [Epump “f ; P\ LkP- 

k=l 


( 7 ) 


The above master equation can be converted into the 
equivalent c-nuinber Fockker-Planck equation in P rep¬ 
resentation, which may be written as a completely equiv¬ 
alent stochastic differential equation [l^ 

+ ( 8 ) 

where a = [a^, Osi, ccii, as 2 , Q:i 2 , a*, a*]^, a*^, 0 * 2 , 0 * 2 ]^, 
and F is the main part of the system’s evolution, in the 
form of = [/,/*]^, where 


/ = 


/ e - 'jpap - 2ga*asiaii - ga*j^as 2 aii - ga*iai 2 asi + galia *2 + gafia *2 \ 

bslCFsl -f gCX^OCj^-^ -f gOii\OCp(y.^2 9^s2^il^p ‘^9 ^P^s2^s1 

bilCFjl -f gCXpCXg^ g(XslOipCyg2 

~7s2Q!s2 + go'siO!pa*i -f ga^^a* 

\ -7i2ai2+5aiiapa*i Tgafittp / 


Matrix B contains the coefficients of the noise terms thus 
it is not included in the stability analysis. To obtain B 
we hrst introduce diffusion matrix 


where d is a matrix 


d = 


—gas20iii 

—gctsl0!i2 

go^h 

V 


-gas2aii -gasiai2 ga^^ gaf^ 

‘^go^p^s2 g^p ^ g^ii^p 

goip 2gotpOii2 gois\ocp 0 

0 gasiap 0 0 

ganUp 0 0 0 


D = 



By dehnition BB"^ = D, then if the diffusion matrix 
is positive semi-dehnite, we can obtain B through fac- 
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toring D, which is the case in our numerical calcula¬ 
tion. rj = where r], 

are real noise terms characterized by {rii{t)) = 0 and 


IV. LINEARIZED QUANTUM-FLUCTUATION 
ANALYSIS 


To solve EqlHl we decompose the system variables 
into their steady-state (classical) values and negligible 
quantum fluctuations around the steady-state values as 
ai = Ai + Sai- Then one can use linearization analy¬ 
sis as a method to calculate the spectra for the output 
cavity modes. The signal photons (sl/s2) share similar 
quantum characteristics with their corresponding idlers 
(il/i2). Because uJk {k = p, si, s2, if, i2) are approxi¬ 
mately equal, in order to simplify the calculation, we 
assume all modes have the same intra- and extra-cavity 
damping rate in the cavity (i.e., 7 fc = 7,7fcc = 7c,7fco = 
7 o, k = p, si, s 2 , zl, i2) in the following part of this paper. 

First we calculate the steady-state solution, so dct/dt 

I 


in Eq. | 8 ]is set to zero. The pump threshold is 

( 9 ) 


When e < eth, the steady-state outputs can be obtained 
as Ap = e /7 and Aj = 0(z = zl,sl,i2,s2). When e > eth, 

the steady-state outputs becomes 


A _ s+y/e'^ + {3'y^)/g 


A^l — -Asl — An — 


Ab = 


isl 

2gApAl 
7 


37 

Y^I^(l - and A ,2 = = 

In the present scheme we only consider 
the situation for the field modes to oscillate above the 
threshold. Now that the 7 is of the order of 10®, Ap = 

— < e /7 is actually much smaller than e, so our 
non-depletion assumption is self-consistent. 

Then we minus the steady-state part from both sides 
of the original equation. The fluctuation part is thus 
linearized as 


A 


5a. = M5a + Brj, 


( 10 ) 


where 5a = [5ap,5asi,5aii,5as2,5ai2, H.cA. M is the 
drift matrix with the steady-state values inserted, given 

by 


M = 


mi m 2 
971,2 


-7 


~gAaAb — 2gApAa 

+ 2-gAaAb 

~gAaAb — 2gApAa -1- 2gAaAb 

- 9^1 

- 9^1 

gAaAh -\- 2gApAa 

- 2gAaAb 
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0 
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0 
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~ ‘2gAaAb 
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0 
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^gApAa 

0 
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0 
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m 2 = 

-gUaUb gal 

‘IgOLpOLb ^Ol-dOLp 

0 






gal 0 

gOiaOtp 

0 

0 






^ g^a gaaap 

0 

0 
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For the linearized quantum-fluctuation analysis to be 
valid the fluctuations must remain small compared to 
the mean values and the eigenvalues of the drift matrix 
M must have no negative real part. If the requirement 
that the real part of the eigenvalues of —M stay pos¬ 
itive is satisfied, the fluctuation eq uations will describe 
an Ornstein-Uhlenbeck process |35l| for which the intra¬ 
cavity spectral correlation matrix is 

S{uj) = (-M -h ia;7)-^£)(-M^ - iwI)-\ (II) 

All the correlations required to study the measurable ex¬ 
tracavity spectra are contained in this intracavity spec¬ 


tral matrix. 

In order to investigate multipartite entanglement, and 
to show that the system under consideration demon¬ 
strates effective five-partite entanglement, we first define 
quadrature operators for each mode as 

Xk = ak + a\, Yk = -i{ak - aj)), (12) 

therefore [Xk,Yk\ = 2i. V{A) = (A^) — (A)^ denotes the 
variance of operator A. 

The output fields is determined by the well-known 
input-output relations [s^. In particular, the spectral 














variances and covariances have the general form 


V. FIVE-PARTITE ENTANGLEMENT 
CRITERIA 
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5£*(w) = 1 + 27,5 x.M (13) 

( 14 ) 


Similar expressions can be derived for the Y quadratures. 

I 


The condition proposed by van Loock and Furu- 
sawa [ 3 ^ , which is a generalization of the conditions for 
bipartite entanglement, is sufficient to demonstrate mul¬ 
tipartite entanglement. We now demonstrate how these 
may be optimized for the verification of genuine hve- 
partite entanglement in this system. Using the quadra¬ 
ture definitions, the five-partite inequalities, which must 
be simultaneously violated, are 


5(1) = ViXp + A,i) + V{-Yp + Ysi + g^Y.i + g,2Ys2 + g^2Y^2) > 4, (15) 

5(2) = V{Xp + Wi) + V{-Yp + 3 , 1^1 + + g,2Ys2 + g^2Y^2) > 4, (16) 

5(3) = U(A,i - A, 2 ) + V{gpYp + F^i + g.^Y.i + y ,2 + g^2Y^2) > 4, (17) 

5(4) = V{Xa - X,i) + V{gpYp + 3 , 1^1 + Fa + g,2Ys2 + Ya) > 4, (18) 


where the guik = p, si, s 2 , fl, i 2 ) are arbitrary real pa¬ 
rameters that are used to optimize the violation of these 
inequalities. It is important to note that in the uncor¬ 
related limit these optimized VLF criterion approach 4. 
Hence, without optimization, some entanglement which 
is present may be missed. 

VI. OUTPUT FLUCTUATION SPECTRA 

From now on, we numerically calculate the values of 
VLF inequalities according to the results obtained above. 
We assume that the resonator is a spherical CaF 2 cavity. 
Note that the theoretical model and analysis are univer¬ 
sal and can be easily applied to other device platforms. 
For CaF 2 , the refractive index is no = 1.43, Kerr coeffi¬ 
cient is 712 = 3.2 X 10“^° m^/W. We assume the CaF 2 
resonator has a radius R of 2.5 mm and light is criti¬ 
cally coupled to the device with a loaded quality factor 
Qo = 3 X 10® (corresponding to a central modal band¬ 
width Awq = 7 p Ri 27r X 64 kHz [ll|), with a pump 
launched at a wavelength of Aq = 1560.5 nm. 

Due to the symmetry relation between signal and idler 
photons, EgfT^ and EqUHlare equivalent. So are EalTTI 
and EqllHl Therefore, we only need to calculate 5(i) and 
^( 3 )- 

To begin with, it is evident from the Eq. [5] through 
m that the intracavity entanglement is completely de¬ 
termined once total damping rate 7, coupling coefficient 
g and pumping power e is fixed, but it is always the vari¬ 
ables outside the cavity that we observe, so the transfer 
also plays a role in the observation, which is mathemat¬ 
ically determined by a ratio, 7 c/ 7 . Thus we first fix the 
first three factors and vary the ratio to see how it affect 
the observed entanglement. In EigUl we plot the mini¬ 
mum of the variances versus the analysis frequency nor¬ 
malized to 7 when 7c takes a portion of 0.34, 0.57, 0.8,1 


of the total damping rate. The blue dashed lines stand 
for 5 ( 1 ), while the green solid ones stand for 5 ( 3 ). 

It can be inferred from the plot that when 7 c = 0 . 347 , 
there is no entanglement between any two of the field 
modes. As we increase the out-coupling coefficients, the 
si and il begin to entangle with the pump photons 
around the center frequency, but it is not until when 
7 c /7 = 0.57 that the s2 and i2 begin to entangle with 
si and il, respectively. Eventually the variance converge 
to EigO (d). Thus we conclude that the entanglement 
among output modes are improved as the 7 c /7 ratio in¬ 
crease, i.e., the entanglement is better when the cavity 
has higher Q therefore lower intracavity loss, and higher 
extracavity coupling coefficient. This can be interpreted 
naturally if we see the coupling as a beam splitter which 
extract squeezed quantum noise to the output [s^, so the 
higher portion the coupling coefficient takes in the to¬ 
tal damping rate, the less consumed entangled pair of 
photons are wasted in the internal loss. Eor that consid¬ 
eration, we will ideally fix 70 = 0 in the following analy¬ 
sis, so that the effect of output transfer is suppressed to 
minimum. 

The pump power, on the other hand, is more subtle. It 
can be derived from our previous analysis that once e/cth 
and w /7 is fixed, the variance is totally determined. In 
other words, the variance Si as a function of w /7 is solely 
determined by the parameter ejeth instead of g, 7 , or e 
independently. This might come counter-intuitive in the 
first sight, as we are indicating that the noise spectrum 
will look just the same except for a scaling in the fre¬ 
quency axis if we the half 712 and 7 simultaneously, as 
long as we keep e/tth constant. However, we should also 
notice that the bandwidth Aw when five-partite entangle¬ 
ment is witnessed will also be halved, indicating a longer 
entangled period in the time domain. When the pump 
is too weak (less than the threshold), there will be no 
sideband, at least not all of them. If the pump grows too 
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FIG. 3. Four variance versus frequency of pump plots when 7 c is 0.34, 0.57, 0.8,1 times as great as 7 (from left to right and 
top to bottom). 7 = 4.02 x 10®s“^, g = 2.21 x 10“"^s“^, e = 1.15eth = 1.97 x lO'^^s”^. The blue dashed curve stand for S(i) 
and S( 2 ), whereas the green solid ones stand for S( 3 ) and S( 4 ). The pump power is fixed at 1.15eth- 
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FIG. 4. The minimum variance as a function of pump power. 


strong, additional sideband that we didn’t take into con¬ 
sideration might pop up, deteriorating the effectiveness of 
our model. We plot the minimal variance throughout the 
noise power spectrum as a function of the pump power 
(normalized by eth) in Fig. 2] and six typical spectrums 
in Fig. [SJ 

It can be inferred from the graphs that both variance 
first descend as the pump power power increases, then 
ascend. S'( 3 ) and 5 ( 4 ) reaches their global minimum at 
e = l.lScth- Considering that they are the short slabs of 
the whole entanglement system, we conclude that 1.15eth 
is the optimal pump power. The other turning point in 
Fig. |3]is around l.leth, when, as we can see in Fig. [S] (b) 
and (c), the variances in the center frequency begin to de¬ 
crease dramatically and become the minimum which was 
once achieved in the side band, as showed in Fig. [5] (a). 


VII. CONCLUSIONS 

In conclusion, we presented a theoretical model for the 
five-partite continuous-variable entanglement among five 
field modes based on cascaded four wave mixing process. 
By solving Fokker-Planck equation in P representation, 
we analyzed the entanglement condition when van Loock 
and Furusawa criteria are violated. We presented the 
design parameters for experimental purpose, and they 
might also be utilized to build integrated compact five- 
partite entanglement generator. We found that the de¬ 
gree of entanglement was totally determined by w/7, 
e/eth and 70 / 7 . This result filled the blank in theory 
for the entanglement analysis of OFCs generated from 
high-Q resonator, therefore would pave the way for fu¬ 
ture optical quantum computation on chip. 
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